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We study the correlations between quark-antiquark pairs in different quantum number channels 
in a deconfined plasma by using an effective model of QCD. Using the three flavour PNJL model, 
the finite temperature spectral functions for different mesonic states are studied at zero and nonzero 
quark chemical potentials. It is found that in the r\ channel resonance structures survive above the 
chiral transition temperature T x , while the kaonic states seem to get washed off just above T x . The 
sensitivity of the structures to the anomaly term are carefully investigated. 
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I. INTRODUCTION 

At high temperatures and densities, strongly interacting matter is expected to undergo a 
transition to new phases where colour degrees of freedom are deconfined. Such conditions of 
high temperatures and densities can be created in the laboratory by the collision of heavy ions 
at large energies. A large amount of data have already been obtained from RHIC. In near future 
results from LHC, for low density and high temperature phase, will be available as well. On the 
other hand, the FAIR experiment at GSI will provide us information about extremely dense 
matter at low temperature. It is thus extremely relevant to study the properties of strongly 
interacting matter at high temperature and density. 

A first principle study of hot and dense strongly interacting matter starting from QCD is not 
easy because the physics is nonperturbative in the temperature and density range of interest. 
At finite temperatures, numerical studies on lattice provide the most reliable results for the 
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physics of the deconfmed phase. The physics of excitations in real time, however, is not directly 
accessible from lattice and one requires an analytical continuation from Euclidean time. Studies 
at finite baryon number density are prohibitively difficult on the lattice, because of the well- 
known sign problem. Some progress has been made in recent years in extracting results for 
small densities [H, [2] . As it stands now, the lattice predicts that the deconfinement and the 
chiral symmetry restoration happen roughly simultaneously Q via a smooth crossover [ij at 
zero quark chemical potential, and the two transitions roughly coincide at about T x ~ 180 — 200 
MeV 31. There is an indication of a critical point at a small chemical potential /is ~ 300 - 500 
MeV, with a first order transition for larger \xb There have also been suggestions that the 
first order line may not exist [2j. 

In view of the limitations in studying finite density and real time problems, it may be useful 
to get an idea of the features of the new phases of matter by studying QCD inspired models. 
Depending on the physics of interest, one can attempt to study a model which has features 
relevant to the physics issue. In particular, for issues of the transition of the hadronic matter to 
a chiral symmetry restored, deconfmed phase the Polyakov loop extended Nambu - Jona-Lasinio 
(PNJL) model has been used. 

In the Nambu - Jona-Lasinio (NJL) model [H, H|, as it is used currently, the interactions 
between quarks are taken into account by suitable four-quark terms which respect the chiral 
symmetry of the QCD lagrangian. Since the gluons are integrated out, with their only effect 
being the four-quark interaction terms, there is no confinement per se in this model. The chiral 
symmetry is, however, spontaneously broken by qq condensation: Ul {Nf) x Un(Nf) — > Uv(Nf). 
As a result, the quarks pick up a constituent mass. Furthermore, the pseudoscalar mesons n, K 
and r\ become the Goldstonc modes. It is also easy to incorporate the V a{X) anomaly by 
introducing a suitable determinant term Q . The NJL model has been widely used to study the 
chiral phase transition in QCD, and also the nature of the excitations at temperatures slightly 
above the transition temperatures. Based on such studies, it was suggested sometime back that 
nontrivial strong correlations between qq pairs persist in the deconfmed phase at moderately 
high temperatures @. 

Due to the lack of confinement, though, the NJL model does miss out on some important 
aspects of the QCD thermal transition. In particular, the chiral transition in QCD is also 
of deconfining nature, as the Polyakov loop, which is the confinement-deconfinement order 
parameter, shows a rapid change. The Polyakov loop extended NJL (PNJL) model [7[ attempts 
to capture this feature of the QCD transition. In this model, the gluon dynamics is described by 
a background temporal gluon field which is coupled to the quarks by the covariant derivative. 
The value of the background field is determined by minimizing the corresponding potential 
[/($) which depends on the traced polyakov loop $. 

The PNJL model has been extensively used to look at the phase structure at high temperature 
and density. Ratti et. al. [|[ have studied the two flavour version of this model in certain amount 
of detail. They have looked at the crossover transition around 200 MeV and also the quark 
number densities at different baryonic chemical potentials. Their results match nicely with the 
lattice, especially the fact that the chiral restoration and the deconfinement transition take 
place almost simultaneously. The speed of sound and the diagonal susceptibilities, calculated 
with the PNJL model, were also found to be in good agreement with lattice results [j| . On the 
other hand the off-diagonal susceptibilities in two flavour PNJL model differ from the lattice 
results Recently Fukushima has studied the phase diagram of the three flavour quark 

matter within the framework of the PNJL model \Tu . 



Since the PNJL model successfully reproduces a large number of quantities calculable directly 
from QCD one may be interested to use this model to study quantities of experimental interest 
which cannot be directly studied from lattice. In particular, a knowledge of the nature of 
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excitations of the plasma is phenomenologically very important. The complete in- medium 
behavior of the mesonic excitations require a study of the spectral functions in these channels. 
As we discussed earlier, it is not easy to calculate them on lattice. A detailed study of such 
spectral functions in the NJL model has been carried out sometime back Q. It is clearly of 
interest to find out how the introduction of the Polyakov loop changes the nature of the qq 
correlations. Hansen et.al. [12] have studied the pion and sigma correlation functions in the 
2-flavour PN JL model. Costa et. al. [l3[ have calculated the masses of the pseudoscalar mesons 
at finite temperature and zero density. It is, of course, interesting and phenomenologically 
important to study a larger set of mesonic correlations for the realistic case of 2+1 flavour 
PNJL model and also at finite density. 

Here we aim to carry out such a study. We look at the spectral functions and pole masses 
of the pseudoscalar channels with different flavour contents. We also carefully investigate the 
effect of the anomaly term in the spectral functions. The spectral functions have also been 
studied in presence of a small chemical potential. 

This paper is organised as follows: in the next section we briefly describe the 2+1 flavour 
PNJL model, as used by us. In Sec. IIIII the basic in-medium calculations and the chiral 
transition have been discussed. The main results of the paper, about mesonic excitations, are 
in Sec. IIV1 A summary of our results and conclusions are available in the last section. 



II. THREE FLAVOUR PNJL MODEL 



In this section, we briefly describe the PNJL model, and specify the parameters of the model 
that we have used. More details can be found in the literature @,B Ell- 
in the NJL model the gluon dynamics is reduced to the chiral point couplings between quarks. 
The PNJL model also introduces the temporal gauge field, since the Polyakov loop, 

* = — J d 3 x jjr Tr Vexp[iJ^dTAi(x,T)] (1) 

is the order parameter for the confinement-deconfinement transition. Here = iAo is the 
temporal component of Euclidean gauge field, j3 — 1/T, and V denotes path ordering. Our 
effective SU(3)f lagrangian is [3] 

£ = E ^/7^V/-E to /^/^/+E^/V'/ + y E [$*v) 2 + (^7 5 a» 2 ] 

f=u,d,s f f o=0,. ..,8 

-g D [dehj> f (l + 75%' + de%(l - Tfe)^/'] - T) 
— Co + Cm + + £ s + £k.mt (2) 

In the above Lagrangian Co is the kinetic term with gauge field interactions, = — iA^S^. 
The gauge coupling is absorbed in the definition of A+ Cm is the mass term which breaks 
the chiral symmetry explicitly. The mass of a particular flavour is denoted by and the 
corresponding field is ipf. The light quark mass is considered to be ~ 5MeV and the strange 
quark mass is considered to be ~ 140 MeV. The term C s is responsible for the four-fermi 
interaction among the quarks with coupling g s . Here we take this coupling to be positive. The 
next term Ckmt which is a six-fermi interaction is invariant under SU(3) L x SU(3) R but breaks 
U{1)a symmetry. This term represents the axial anomaly of QCD Q. Here 'det' stands for 
the determinant with respect to the flavour indices, and the anomalous coupling is represented 
as go- This anomaly term is responsible for the flavour mixing of rjo and r]g mesons in the 
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pseudoscalar channel, giving rise to the r] and rj mesons. The potential is expressed in 

terms of the Polyakov loop $ and its conjugate, $, as @. 

W($,$,T) b 2 (T)- b 3 , o -, N 6 4/ - x2 

-^ = -Y w "j( $ + $ ) + j( w ) (3) 

where 



6 2 (T) = ao + ai (^) + a2 (^) 2 + a3 (^) 



(4) 



b 3 and 64 being constants and To, which is a parameter here, is chosen to be 190 MeV. 

Both $ and $ are treated as classical field variables. When the quark number density is 
zero one has <J> = $. This quantity can be considered as the order parameter for the phase 
transition. Furthermore £/[$,$, T] has a Z(3) center symmetry which encompasses the phase 
transition in QCD. At low temperature U has a single minimum at $ = 0. At high temperature 
U has three degenerate minima at <£> = 1, e ±2i7r / 3 . 

To study the chiral transition, we study the system of eqn. ([2]) in the Mean Field Approxima- 
tion (MFA) to get the field equations for $, $, a. Due to the breaking of SU(3) L x SU(3)r 
symmetry to SU(3)v the quark condensate acquires a non-zero vacuum expectation value given 

by 

< $fl>f >= ~iN c £t y ^ x+ (Tr S f (x - y)) (5) 

where trace is over colour and spin states. Also, eight goldstone bosons appear for Nf = 3 
model. The self-consistent gap equation for the constituent masses are 

M u = m u - 4gs<T. u + 2g D a d a s 
M d = rrid- 'igso-d + 2g D a s a u 

M s = m s - 4g s cr s + 2g D a u a d (6) 

Here Of =< ifcfipf > denotes the chiral condensate of a quark with flavour /, where / = u,d, s. 
The expression for aj at T = and /1 = can be written as [3] 



\/p 2 + M f < 

A being the three momentum cut off. 

The parameters of the NJL part of the Lagrangian are fixed at T = fJ, = 0. In the present 
work we have used the parameter set obtained in ref. @, 

A = 631.4MeV, m„ = m d = 5.5MeV, m s = 135.7MeV, g s A 2 = 3.67, g D A 5 = 9.29 

This parameter set reproduces the physical values of the pion mass, = 138 MeV, pion decay 
constant, / w = 93 MeV, and masses of the K and 77' mesons, ttlk = 495.7 MeV and m,,/ = 957.5 
MeV. The parameters of [/($) are[l [lj], 

a = 6.75, at = -1.95, a 2 = 2.625, a 3 = -7.44, b 3 = 0.75, 64 = 7.5, T = 190MeV 



so as to reproduce the lattice measurements of the Polyakov loop expectation value. 
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III. CHIRAL TRANSITION IN PNJL MODEL 



In the mean field approximation, the thermodynamic potential of the PNJL model can be 
written as [l4| 

n($,*,M,T, fi )=U[$,*,T\ + 2g s J2 <?f 2 -±9D<Ju<J d <J s -TY, I ^Trln^'"^" /,) 



f=u,d,s r 

where S~ x is the inverse quark propagator, 

S- 1 (p ,p)=^ {p° + fji-iA 4 ,)-j-p-M 



(2tt)S 



T 



(8) 



(9) 



and oj n = irT(2n + 1) are the Matsubara frequencies of fermions. Using the identity Trln (X) 
lndet (X), and performing the colour trace, we get 



f—u,d,s 



d 3 p 



V Jo (2tt) j 



In 



l + 3($ + $e ^ ) e r 



In 



l + 3($ + $e- 



(g p/ +M) (g p/ +M) 



(10) 



where E^/ = y P 2 + -My is the single quasiparticle energy. In the above integrals, the vacuum 
integral has a cutoff A whereas the medium dependent integrals have been extended to infinity. 

To look at thermal behavior of different obscrvables, we need to study the variation with 
temperature of the minimum of the thermodynamic potential. In fig. [T] we show the thermal 
behavior of the chiral condensates for both the light and the strange quarks, in both NJL and 
PNJL models. The NJL results can be obtained from eq. ([TO]) by setting [/[$, <&, T] = 0, $ = 1 
and $ = 1. Fig. [T] shows a rapid drop in (uu) in both the models in the temperature range 
150-200 MeV, without any singularity, indicating a chiral symmetry restoring crossover. This 
is consistent with what lattice studies find for QCD [J]. From the points of inflexion of (uu)(T) 
one estimates T x - 170 MeV for NJL and 190 MeV for PNJL models. Figure Q] also shows 
the behavior of the Polyakov loop, $ [lj|. The figure shows that deconfinement temperature, 
Td ~ T x , as expected from lattice studies of QCD Q. 

At large quark densities, the chiral symmetry restoring transition in PNJL model becomes 
first order. Figure [T] also shows the behavior of the chiral condensate and the Polyakov loop in 
the presence of a chemical potential /i 9 = 320 MeV, which is just above the critical point for 
this model. The transition temperature at this /i 9 is ~ 85 MeV. Note the sharp fall of (uu}, 
indicating a first order transition. All the quantities shown in the figure show a jump at the 
transition point, consistent with a discontinuous transition, though the drop in (ss) and the 
rise in $ are only a few %, and both these quantities have rather complicated temperature 
dependence above the transition point. In the presence of a chemical potential, there is no 
charge conjugation symmetry so $ 7^ as seen in fig. [TJ 

One problem of the mean- field study of PNJL model is the absence of color neutrality at finite 
density [13, E|- This arises because of the fact that the Polyakov loop, which couples to the 
dynamical quarks, serves as an external coloured field, leading effectively to separate chemical 
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FIG. 1: (Color online) (a) (uu)(T)/(uu)(T = 0) (red, solid) and <ss)(T)/ {ss)(T = 0) (green, dashed), 
as a function of T/T x at /j, q = for the NJL model, (b) The same for the PNJL model; also shown 
is the Polyakov loop expectation value, <3> (blue, dash-dotted), (c) Same as (b) but for ji q — 320MeV. 
Here $ (blue, dash-dotted) and $ (pink, dotted) are different. 



potentials for quarks of different color. As a remedy to this problem, it has been suggested to 
introduce a colour chemical potential which will enforce color neutrality pjj. While this affects 
the individual number densities of the quarks of different colors, the total number densities of 
quarks and antiquarks, n q and nq, do not change significantly We will discuss in Sec. IIVBI 
the implication of this for our study. 

In fig'H] we have shown the variation of M u and M a with temperature at different chemical 
potentials. Here M u and M s are the constituent quark masses, eq.©. At ^ q = 0, The change 
in the mass is smooth. But the change becomes sharper on introducing a chemical potential. 
While the introduction of the Polyakov loop changes the constituent quark masses considerably 
around the phase boundary, at very high temperature both the models give similar results, 
since there $ — > 1. A similar variation of the strange quark mass is obtained in the two models. 

The phase diagram we obtain is shown in fig. [3] The critical point in our calculation was 
found to be at (/j, c ,T c ) — (314 MeV, 92 MeV). While lattice estimates of the critical chemical 
potential fi c vary considerably between different groups, they tend to be much lower than the 
PNJL value, < 150 MeV. This probably indicates that the naive PNJL (or NJL) model may 
not be a good model for the QCD transition at large fj, q ~ 300 MeV, since nuclear excitations 
become important in this regime. 

Another point to mention here is that the L0t(l) anomaly term is introduced explicitly in 
the NJL and PNJL models by a six fermion interaction term. In QCD, of course, this term 
is introduced due to instantons, which are topological objects in the gauge sector [ijj]. It 
is not known how well the anomalous Ua(1) symmetry is restored in the deconfined phase. 
In Ref. [6j it was speculated that the Ua(X) symmetry is restored well below T x . While 
there is no evidence of that from lattice studies, it is an open question whether the Ua{^) 
symmetry is restored (or at least, the effect of the anomaly term substantially reduced) at 
above T x . If we include an explicit temperature dependence of the coefficient of the anomaly 
term, go(T) = .9d(0) exp(— T/T x {mu)) 2 , where T x (^jl) is the transition temperature at quark 
chemical potential /i, then the transition becomes considerably softer and the first order line 
completely disappears, while the transition temperature moves to a somewhat lower value (see 
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FIG. 2: (Color online) (a) Variation of M u (a) with temperature, at /i = and 320MeV for both NJL 
and PNJL model. Solid (red) and dashed (green) lines are for NJL model with fi q — and 320 MeV, 
respectively, while dashed-dot (blue) lines and dotted (pink) lines show the corresponding figures for 
PNJL model, (b) Same as (a) but for M a . 
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FIG. 3: (Color online) Phase diagram in the PNJL model. The solid line at large fi q ends in a critical 
point at (fi c ,T c ) — (314 MeV, 92 MeV). At smaller \i q one gets a crossover, denoted by the dotted 
line. Putting in an explicit temperature dependence in the anomaly term go{T) = go(0) exp(—T/T x ) 2 
(denoted PNJL1; see text) washes off the first order line completely. 



fig. [3]). We use PNJL1 to denote the model with such a temperature dependence of go in eq. 
© ■ One might be tempted to make a conclusion that the existence of critical point is strongly 
related to the strength of the instanton- induced anomaly term at the transition point [2fj| . 

Of course, PNJL1 implies a substantial reduction of the anomaly term already before T x . A 
more realistic temperature dependence of go could be a sharp drop around T x . To investigate 
the effect of such a temperature dependence, we also consider eqn. © with <?r>(T) = go{0) (1 — 
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tanh 5x)/2, where x = {T -T x {p)) /T x {p). We call this model PNJL2 [U. PNJL2 has a phase 
diagram similar to PNJL i.e. the critical point has the same position at (/x c ,T c ) = (314 MeV, 
92 MeV). However, as we will see later, the nature of the mesonic excitations around T x can 
be substantially different between the two. 



IV. MESONIC EXCITATIONS AND 17a (1) ANOMALY 

In order to understand the properties of the medium beyond the bulk thermodynamic prop- 
erties, we need to look at the low-lying excitations. The simplest gauge invariant excitations 
in QCD are the mesonic excitations. Masses and decay widths of the mesonic resonances are 
calculated from the correlations of grq-type operators in QCD vacuum. In the medium, we are 
interested in spectral changes due to interactions with the medium; in particular, we would be 
interested in possible resonance-like structures in the deconfined phase. We need therefore to 
look at the spectral functions. 

Here we study spectral functions of pseudoscalar and scalar mesonic states. They are of 
interest because of their close connection with the chiral symmetry breaking and its restoration. 
The temperature dependence of the spectral function has been studied in the Nambu-Jona- 
Lasinio model in ref. @. There it was found that narrow structures persist in the symmetry 
restored phase, at moderately high temperatures. Here we intend to see how the coupling with 
the Polyakov loop affects these results. Also we study the sensitivity of these structures on 
the coefficients of the anomaly term. We also extend the studies to nonzero quark chemical 
potentials. 



A. Formalism 

The spectral function um (lj, k) for a given mesonic channel M in a system at temperature T 
can be defined through the Fourier transform of the real time two point functions D > and D < 

a M {u,k) = ±-{D>(k ,k)-D>(k ,k)) (11) 



D>(x ,x) = (J M (x o ,x)J M (0,0))c (12) 
D^(x ,x) = (Jm(0,0)Jm(xo,x)) c 

Here Jm is the suitable hadronic operator and the subscript c denotes the connected part of 
the correlation function. Using the Kubo-Martin-Schwinger(KMS) condition [22] 

D>(x ,x) = D < (x +i/T,x) (13) 
one can connect itm(w, k) to the retarded correlation function, 
(j M {u,k) = 2ir 1mDff(ko,k) 



D«(k ,k) = I j^e«-*0(x o )([J M (x o ,%,JM{O,d\)c (I 4 ) 
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Inserting a complete set of states in eq. ([TTT) and using eq. (fT3|) . one gets the expansion 

CTM ( w ,fe) = ^^(e-^/ r ±e-^/ T )Kr 1 |JM(0)|m}| 2 5 4 (^-^ + C) (15) 



in , n 



where Z is the partition function, and q n refers to the four- momenta of the state |n). 
A stable mesonic state contributes a 5 function-like peak to the spectral function: 

a M {u,k) = \(0\J M \M)\ 2 e(k Q )S(k 2 - m 2 M ) (16) 

where mu is the mass of the state. For an unstable particle one gets a peak with a finite width, 
e.g., the Breit-Wigner form for states with narrow decay width. We want to study how the 
spectral function changes due to collisions with the thermal medium. 

At the level of approximation we are working, the collective excitations, the fluctuation of 
the mean field around the vacuum, can be handled within the Random Phase Approximation 
(RPA) [23| . In this approximation, which is equivalent to summing over the ring diagrams, the 
retarded correlation function is given by 

Here G is the suitable coupling constant and IlM(fc 2 ) is the one-loop polarisation function for 
the mesonic channel under consideration, 

HM(fc 2 ) = J ■^ i Tr{T M S(p+^)T M (p-^)}. (18) 

S(p)is the quark propagator and Tm is the effective vertex factor. Eqn. (|18j) has been evaluated 
in the literature; for the NJL model, the effective formulae were calculated in ref. @ . The whole 
effect of the introduction of the background Polyakov loop can be absorbed into a modification 
of the Fermi-Dirac distribution functions 

(l> + 2^e~ f}{ - Ep+ ^ ) e r^ E p+^ + e -3P( E p+v) 
/<£> (Bp) = 1 + + $e _ /3(Bp+M )- )e _^ (Bp+Al) + e - 3 /J(is P + A1 ) ( 19 ) 

($ + 2$e~ /3(Bp ~' j) )e~^ ( - B p~^ + e ~ 3 ' 3( - E p~^ 

/$ ( E p) = ! + 3 ($ _|_ $ e -/3(£ p - M )) e -/J(£p-M) + e -3P(E p -v) ( 20 ) 

For NJL model, $ = 1 and the distribution functions become the usual Fermi-Dirac distribution 
function. 

From the correlation function one can calculate the spectral function using cq. (fT4|) . The 
position of the pole controls the decay of the correlator and is called the pole mass. It can be 
obtained by solving 

1 - 2Gn M = 0. (21) 

The coupling constants G for the different flavour combinations are given in table 1. The 
cases of r\ and ?/ need a special mention. These mesons arise from the mixing of the flavour 
singlet and octet states, i.e. rjo and 77s mesons. The mixing arises because of an interplay 
between the anomaly term and the nondegeneracy between m s and m Mj( j. The masses of 77 and 
rj' may be obtained from the roots of the equation 

det[l — 2IIG] = det \ (^ncroC* 00 + 2n 80 G8o) — (2n oGo8 + 2iI 8G88) | — q (22) 

\ — (2IIooG80 + 2Il88G8o) 1 — (2il8oGo8 + 2II88G88) / 



10 



G„ = gs + gD-y 


G K ± 


= gs +9dP 


G K ° = gs + gooL 


Goo = gs - 


- Z(a + /3+~/)g D 


G33 = Gtt 


Gss = 98 - 


i( 7 -2a -2)8) ff D 


G&=—$ s (a-0)gn 


^ss — 











TABLE I: Pseudoscalar coupling strengths in SU(3) NJL model, where a =< uii >, /3 =< dd >, 

7 =< ss > 

where 

n 00 = |[2iC(fc) + n s p s (fe)] 
n 08 = (^)[nl(fc)-nf s (fc)] (23) 
n 88 = (~)[iC (fc) + 2nf s (fc)] 

The polarisation functions are obtained from eq. (|18l) . On the other hand pole masses of 770 and 
r/ s can be obtained directly from n o, n 88 and eqn. (l2"Tj) . 

B. Results at finite temperatures 

The pion and sigma channels provide the most interesting observables with mesonic states, 
since they are directly associated with the chiral symmetry. In the chiral symmetry restored 
phase, pion and sigma should be degenerate. Below T x , on the other hand, pion is the goldstone 
mode of the spontaneously broken chiral symmetry and is therefore much lighter than the sigma. 

In fig. U]we show the temperature dependence of the pole masses of the pion and the sigma in 
the NJL and the PNJL model. We see that close to T x there is substantial difference between 
the two models. While the two states come close to each other, they are not degenerate just 
above T x . As mentioned in Sec. IIII1 results of the two models converge at higher temperatures. 
This, combined with the fact that T x is lower for NJL model, leads to a smoother change in 
NJL model above T x than in PNJL. The nondegeneracy between the pion and the sigma can 
also be seen in the fig. 4(b) , which show the spectral functions in the pion and sigma channels 
at a temperature of 1.1 T x . 

There can be two sources for the nondegeneracy: the finite (though small) current quark 
masses, and the remnant /7a(1) symmetry. To investigate the effect of the Ua(^) symmetry, 
we look at PNJL2, which introduces a cutoff in <?£>, as explained in sec. IIIII As can be seen in 
fig. 21 this significantly alters the spectral functions, and removes most of the nondegeneracy 
of the 7T and a above T x . Also, while there is a clear structure in both of these channels above 
T x , the structure is considerably broader in the absence of the anomaly term. When studying 
remnant bound state structures from QCD-inspired models, it therefore is important to keep in 
mind the effect of the anomaly term, since unlike in QCD, the anomaly terms are put in here 
by hand and therefore their contributions across phase transition may be different from that in 
QCD. 

In figure [S] we have plotted the masses of 77, 77', rf and 77 s mesons respectively. We see 
features similar to those in fig. 2) significant differences around the transition region due to the 
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FIG. 4: (Color online) (a) Pion and sigma mass versus temperature, at [i = MeV for both NJL (red, 
solid and blue, dash-dotted) and PNJL (green, dashed and pink, dotted) model, (b) Spectral functions 
for the pion and sigma in the PNJL model, just above T x . Also shown is the effect of suppression of 
the anomaly term at T x (PNJL2; see Sec. IIII|) on the spectral functions . 



introduction of the Polyakov loop. It is also seen that very close to the transition region, the 
masses of the singlet and octet states start approaching each other, indicating that the effect 
of the anomaly term starts getting reduced. However, without an explicit cutoff like PNJL2, 
the nondegeneracy between these two states is clearly visible at 1.1 T x . 
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FIG. 5: (Color online) Masses of the different r? states in NJL(green, dashed) and PNJL (red, solid) 
model. 



The spectral functions of the rjo and rjs states, calculated in the PNJL model, are shown in 
fig. El Here we also investigate the effect of the anomaly term by comparing the calculations in 
PNJL2 model. We find here that the effect of the anomaly term is considerable in 770 but not 
in 77s . This can be understood by looking at the two relevant couplings, Gs8 and Goo, in table 
1. Because of the different numerical factors, and also the different combinations of the quark 
condensates, the anomalous coupling of 770 is substantially larger than that of 77s. Interestingly, 
the broadening due to anomaly suppression is much less even for 770 than was seen for the pion. 
In these channels, a clear (though broad) resonance structure is seen just above T x . 
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FIG. 6: (Color online) (a) Spectral functions of 770,8 at temperatures of 0.9 T x and T x . Note the 6 
function peak for r/s below T x . (b) The same above T x . Also shown are the effects of cutting off the 
anomaly term above T x at T x (PNJL2; see Sec. 



The behavior of the pole mass for the kaon state, fig. [7j is qualitatively similar to the pion. 
However, the spectral function above T x , shown in the same figure, shows a very different 
behavior. The sharp peak at T x is completely washed off already by 1.1 T x : there is no 
indication of a strong correlation between quarks in this channel. 
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FIG. 7: (Color online) (a) Pole mass of the kaon state, in PNJL (solid, red) and NJL (dashed, green) 
models, at different temperatures, (b) Spectral function, from the PNJL model. The sharp peak at T x 
seems to be completely washed off already by 1.1 T x . 



C. Results at finite densities 



One of the advantages of model studies is that one can study regions in parameter space 
which are not easy to study directly from QCD in a controlled way. The agreement of PNJL 
model results with QCD at finite temperatures encourages us to investigate effects of nonzero 
baryon densities, by introducing a quark chemical potential. As mentioned in Sec. IIII1 however, 
we do not expect the PNJL model to be a good model for QCD transition at large fi q . We 
therefore investigate the mesonic excitations in this model at /Lt q = 150 MeV. This is around 
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the region where lattice studies suggest the beginning of the first order line, while the PNJL 
model still finds a crossover transition at these densities. From fig. [3] we also get T x ~ 170 
MeV at fi q = 150 MeV, which is not too far from what one gets from lattice. 

In fig. [H we show the pion and sigma spectral functions above T x for (j, q = 150 MeV. We see 
that the spectral functions have a sharp structure even substantially above T x , as in the case 
for n q = 0. The nondegeneracy of the two channels and the effect of the anomaly terms are 
also similar to the fi q = case. In the figures 8(b) and 8(c) we show the variation of the pion 
pole mass with for two different temperatures. Interestingly, the quantitative difference in 
the masses calculated in PNJL and NJL models increases with density. Of course, as we have 
discussed above, we do not expect these models to be a good approximation to QCD at values 
of n q > 150 MeV. 
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FIG. 8: (Color online) (a) Spectral functions of pion near the chiral transition point, for /i 9 = 150 MeV. 
Also shown is the effect of suppression of anomaly term at T x (PNJL2; see Sec. IIII|) . (b) Pion mass 
versus fj,, at T = lOOMeV. Solid red line is for PNJL model and dashed green line is for NJL model, 
(c) Same as (b) but at T = 150 MeV. 

The spectral functions of t]q and r]$ are shown in fig. [9l for /j, q — 150 MeV. This figure has 
some interesting features. First, the interplay of the temperature and chemical potential gives 
rise to a nonmonotonous behavior in the peak of the rjs spectral function in the deconfined 
phase. The peak position drops sharply at T x , and then rises again. Second, we note that the 
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7/8 peak, which is below the rjo peak below or at T x , becomes very broad already at 1.05 T x , 
with the peak position slightly above that of the i]q peak. If this feature is true for full QCD 
and is not an artefact of the model, it will have interesting phenomenological implications. 
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FIG. 9: (Color online) (a) The spectral functions of the r/o channel around the chiral transition tem- 
peratures, for fi q = 150 MeV. (b) Same as (a) but for rjs- 



Let us now discuss the sensitivity of our results on enforcement of colour neutrality (Sec. 
IIIip . Our calculations at finite density have been carried out at fi = 150 MeV and T = 1.1T C . 
At such T and /i, if we want to enforce color neutrality by the introduction of a color chemical 



potential, we need a color chemical potential + /UgTs (in the notation of Ref. [17|) with 
a small jUs(~ 3QMeV), while /X3 is negligible. Here T3 and T$ are Gell-mann matrices in 
color space. While this changes the individual colour densities, the total quark number density 
(n q — riq) and the total scalar density (n q + rig) do not have any appreciable change. Similar 
results about the total quark density has also been reported in ref. [18[ . In fact the change in the 
scalar density is less than 5% under such conditions. Here we have studied spectral functions 
of colorless, meson-like objects, which are sensitive to the scalar density. Hence we do not 
expect our results at finite chemical potential to be affected appreciably by the issue of colour 
neutrality. 



V. CONCLUSION 



To conclude, we have studied the meson-like excitations in the deconfined plasma at both zero 
and nonzero quark chemical potential, using the Polyakov loop extended Nambu - Jona-Lasinio 
(PNJL) model. 

The phase diagram of the PNJL model shows a first order line of transitions at large fx q , 
ending in a critical point at (/i c ,T c ) ~ (314 MeV, 92 MeV). We note, however, that the phase 
diagram can be quite sensitive to the temperature dependence of the coefficient of the anomaly 
term. In particular, a suppression of the anomaly term at large temperatures may lead to a 
complete wash off of the first order line. 

The in-medium behavior of the mesonic correlations was studied by looking at the spectral 
functions. In the pion and sigma channels, reasonably sharp structure was found still at 1.1 
T x ; also the two states are not degenerate at these temperatures. While both these features 
can be phenomenologically interesting, it was also found that these features are quite sensitive 
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to the anomaly term. In particular, suppressing the anomaly term above T x leads to much 
broader and degenerate structures in these two channels. Since the structure of the ir resonance 
above T x is of importance for RHIC phenomenology, this sensitivity implores one to look for 
the restoration or otherwise of the U(A)i anomaly at T x by looking at suitable observables 
directly on lattice. 

The 778 and the 770, on the other hand, show much less sensitivity to the anomaly term. Both 
these channels show a clear resonance at 1.1 T x . rjg shows very little sensitivity to a suppression 
of the anomaly term, 770 is more sensitive, but the resonance-like structure remains even when 
the anomaly term is suppressed. In the PNJL model, the kaon channel does not show any sharp 
structure above T x . All these behaviors will be of interest for RHIC phenomenology. 

We also looked at the mesonic spectral functions just above T x for a nonzero fi q . At /j g 
= 150MeV, which should be in the regime covered by the RHIC energy scan, PNJL shows a 
crossover transition. The features of the spectral functions in the pion and sigma channels are 
found to be quite similar to those at [i q = 0. A more interesting behavior was found in the r\ 
channels. The interplay of the finite fi q and finite T lead to a nonmonotonous behavior of the 
peak in the rjs channel. Also it was found that the relative positions of the rjo and 77s peaks 
get reversed above T x : the 770 peak comes at a lower u than the 77s peak at 1.05 T x . It will be 
interesting to explore the phenomenological implications of these behaviours. 
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